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Abstract 

By means of synchrotron X-ray diffraction, we studied the effect of high pressure, P, up to 13 GPa 
on the room temperature crystal structure of superconducting CaCg- In this P range, no change of 
the pristine space group symmetry, R3m, is found. However, at 9 GPa, i.e. close to the critical value 
at which a large T c reduction was reported recently, we observed a compressibility jump concomitant 
to a large broadening of Bragg peaks. The reversibility of both effects upon depressurization and 
symmetry arguments give evidence of an order-disorder phase transition of second order, presumably 
associated with the Ca sublattice, which provides a full account for the above T c reduction. 
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The observation of superconductivity with remarkably high critical temperature T c =11.5 
K in CaCfi renewed the interest in graphite and related compounds as a model system 

with a flexible layered structure and tunable electronic and transport properties. Penetration 
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depth [3j , tunneling [4| and specific heat [5[ experiments and Density Functional Theory (DFT) 
calculations consistently point at a conventional s-wave BCS scenario with a sizable 

electron-phonon coupling constant, A ~ 0.8, that mainly involves two phonon modes with in- 



plane Ca and out-of-plane C polarizations. According to further studies 
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;he whole series 

of AC$ compounds hitherto synthesized containing alkaline earths A=Sr,Ca,Ba 

BflQ, the 

largest T c found in CaC6 arises from a stronger electron-phonon coupling between the interlayer 
band and the above phonons. In addition, the above studies suggest a scaling of T c with the 
inverse of the interlayer distance that can be tuned either by pressure or by the A cation size. 

Consistently with this scenario, T c is found to^ increase with pressure in all ACq compounds, 
i.e. CaC 6 



]J, Q, Q, SrC 6 |s] and YbC 6 [151 ] . CaC 6 exhibits the largest increase up to 



T c =15.1 K, the highest value for graphite-related compounds, at 7.5 GPa 12J. At P ~ 8 GPa, 
a sudden T c drop down to ?^5 K concomitant to an anomaly of the aforementioned Ca and 
C phonons suggests a picture of pressure-induced structural instability limiting T c 12[. This 
picture is supported by the theoretical prediction of a vanishing mode frequency for the in-plane 
Ca phonon under high pressure [s[ [l^. A similar instability seems to occur also in YbC6 at 



P = 2 GPa [15[. Whilst the above DFT calculations provide a consistent description of the 
electronic structure, controversial predictions were made as to the P-dependence of the phonon 



modes and of T r 



0, y, Q 



17| . For instance, the P-dependent A values are systematically 



underestimated considering the large increase of T c found experimentally. Also, the structural 
instability responsible for the drastic T c reduction at 8 GPa and its effect on the relevant BCS 
parameters remain to be determined. A possible reason of these difficulties is the inaccurate 
calculation of the high-pressure structure that remains to be determined experimentally. 

In order to address this point, in this work we carried out a detailed study of the room 
temperature crystal structure of CaCg at high pressures up to 13 GPa. At 9 GPa, i.e. close 
to the borderline between the high- and low-T c phases, we found a jump of the compressibility 
concomitant to a large Bragg peak broadening. We show that these results put into evidence 
an order-disorder phase transition of second order that accounts for the above T c reduction 
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and structural instability at 8 GPa. To our knowledge, this is the first case of superconductor- 
superconductor phase transition driven by an order-disorder structural transition. 

The experiment was carried out at the high-pressure diffraction beamline ID9 of the ESRF. 
A high-quality mm size bulk sample of CaCgwas prepared from a platelet of highly ori- 
ented pyrolithic graphite, as described elsewhere UJ. A few ~50-100 /im pieces were mounted 
in two diamond anvil cells (DAC) in the opposite anvil configuration. The two DACs were 
subsequently charged with high purity helium (sample n.l) and argon (sample n.2), used as 
pressure-transmitting media. Due to the reactivity of CaCg, the procedure was carried out in a 
high-purity dry box, where the cells were sealed under an initial gas pressure of about 0.1 GPa. 

The diffractograms were taken in the Debye-Scherrer geometry suitable for the above DAC 
configuration using a wavelength A = 0.4133 A. The diffracted signal was collected within a 
maximum diffraction angle 2d max = 25° by an image plate. Owing to the small diameter of the 
stainless steel gasket (150 /im), the beam spot size was reduced down to 40 fim for P < 10 GPa. 
At higher P, this size was further reduced to 20 /im due to the shrinkage of the cell. For each 
diffractogram, the P value was determined from the emission line of a ruby crystal. To maximize 
the number of reflections detected, the DAC was spindled within an angle of ±3° with respect 
to the incident beam direction. This is especially important in our case, as the layered structure 
of CaCg is expected to induce a preferential orientation of the graphene layers perpendicular to 
the incident beam. For both samples, no significant amount of secondary phases was detected 
in the sample regions probed by the X-ray beam. This enabled us to reliably refine the unit 
cell up to 13 GPa. Sample n.l gave a larger signal-to-background ratio, so the analysis below 
mainly refers to this sample. More than hundred diffractograms were taken on sample n.l by 
progressively increasing P from 0.1 up to 13 GPa and upon depressurization. Representative 
diffractograms are shown in Fig. 1. All of the main Bragg peaks are found to match those 



of the rhombohedric R3m phase The indexation reported in 



equivalent hexagonal unit cell with a=4.330 A and c=13.572 A This cell choice enables 



'ig. 1 corresponds to the 



to better visualize the crystal structure as an alternate stacking of graphene and Ca layers (see 
Fig. 1). The additional features appearing above 9 GPa are discussed below. 

The indexation of Fig. 1 shows that, despite the expected preferential orientation of the 
sample, reflections with various combinations of h, k and I indices are present. This might be 
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FIG. 1: (Color online). Representative diffractograms of sample n.l at different pressures during and 
after depressurization. Intensities are normalized to the integrated intensity of the strongest (113) 
peak. Bragg indices refer to the hexagonal unit cell in the inset. Diamonds mark unindexed peaks 
and additional broad features appearing above 9 GPa (see text). 

due to the folding of the CaC6 flakes during sample manipulation. About thirty diffractograms 
in the 0.1-13 GPa range were selected for the unit cell refinement. For the P < 10 GPa 
data, we typically used as many as 12-15 Bragg reflections. This number reduces to 8 at 
higher pressures because of the weaker signal and the appearance of additional broad features 
to be discussed below. The result of the refinement and the resulting equation of state are 
shown in Fig. 2. Thanks to the presence of strong peaks, e.g. (300), (113), (003) and (006), 
corresponding to independent Bragg planes, the accuracy obtained for the a and c parameters 
is as good as ±0.001 — 0.002 and ±0.002 — 0.004 A, respectively. Good refinements were 
obtained also using the data of samples n. 1 and n. 2 upon depressurization and pressurization, 
respectively. By comparing these three data sets, one notes that the data are reproducible 
within the experimental error and no hysteresis is found after depressurization. 

Fig. 2 indicates that both a and c parameters decrease monotonically and roughly lin- 
early with P down to 9 GPa. In this P-range, a linear regression yields the following val- 
ues of isothermal compressibility coefficients: da/ dP '=-0.0038 A/GPa; dc/ dP=-0.081 A/GPa. 
Hence, the c-axis compressibility is about twenty times larger than the in-plane one, as 
expected considering the layered structure. The corresponding volume compressibility is 
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FIG. 2: (Color online), (a) Experimental (symbols) and calculated (solid lines) pressure-dependence 
of the hexagonal a and c unit cell parameters, (b) Pressure-dependence of the unit cell volume 
corresponding to the data (symbols) and calculations (solid line) of panel (a). The T c vs. P data from 



ref. 



121 ] are also reported to show the coincidence of the midpoint of the T c drop with the kink of the 



volume compressibility at 9 GPa marked by vertical arrows. The broken line is a guide to the eye. 

k = (1 /VgQp a )dV/ dP = —0.0082 GPa -1 . The above values and the large anisotropy are consis- 
tent with previous pressure studies on graphite 18j and graphite intercalated compounds with 
same composition and similar crystal structure, such as LiC6 [jj]]. In order to study in detail 
the P-induced structural changes, the relaxed structure was calculated as a function of pres- 

n 

sure up to 12 GPa by means of the ESPRESSO implementation I2CJ of the DFT, as described 

on 

elsewhere p, |8|]. The results of the calculations are compared with the experimental data in 
Fig. 2. For P < 9 GPa, it is noted a good agreement for both a and c parameters, except a 
slight overestimation of the former. 

Most notable feature of the experimental equation of state of Fig. 2 is a kink at 9 GPa 



5 



corresponding to a jump of the isothermal compressibility, An, the signature of a second order 
phase transition. The kink leads to a discrepancy between experimental and calculated data in 
the P > 9 GPa range, for, in this range, we performed only volume optimization calculations 
which do not take into account possible structural instabilities. A linear fit of the experimental 
data in the P > 9 GPa region yields k = (l/Vgcp^dV/dP = —0.0215 GPa" 1 , i.e. almost 3 
times larger than in the P < 9 GPa region. In order to interpret this result, we notice that 
no change of space group symmetry is seen in the diffractograms of Fig. 1. Thus, the jump 
of k at 9 GPa suggests a structural change allowed by symmetry leading to a discontinuity 
of the electronic structure. It turns out that the only free parameter of the R3m structure 
is the in-plane x coordinate of the C atoms in the hexagonal cell 2l|. Thanks to this simple 
symmetry argument and considering the rigidity of the C honeycomb layers within the plane, 
one concludes that the only atomic displacement allowed by symmetry is the off-centering of the 
Ca atoms with respect to the C layers. At ambient pressure, the Ca atoms are centered with 
respect to these layers 8], as expected from a simple consideration of energy minimization. 
This implies x = 1/3, whilst the Ca position is (0,0,0). However, an off-centered in-plane 
displacement of the Ca atoms requires a modest amount of energy due to the soft phonon mode 
associated with this displacement j3] . This scenario is indeed supported by a recent DFT study 
by Csanyi et al. on the structural stability of CaC6 at high pressure 22|. 

According to the above picture, the large broadening of all peaks concomitant to the kink at 
9 GPa (see Fig. 3) shows that the in-plane displacements of the Ca atoms are highly disordered. 
The broadening turns out to be intrinsic, as it is completely reversible upon depressurization. 
We conclude that, at 9 GPa, an order-disorder transition of the Ca sublattice occurs. A simple 
empirical account for this transition is as follows. The P-induced compression of the unit cell 
along the c-axis leads to the reduction of the interlayer volume available for the Ca atoms. 
Because of the relatively large size of the Ca ions, this volume reduction further reduces the 
limited stability of the potential minima associated with the symmetric positions. Beyond a 
certain point, any position between adjacent graphene layers becomes energetically equivalent. 
This disordering mechanism is consistent with the previous observation of anomalous hardening 



of the out-of-plane C mode and of 



3ad metallic properties with high residual resistivities and 



flat resistivity curves above 8 GPa 121 ] . 
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FIG. 3: (Color online), (a) Pressure-dependence of the full width at half maximum (FWHM) of the 
(113) (circles) and (300) (squares) peaks. Open symbols refer to data after depressurization. Error 
bars are smaller than symbols when not visible, (b) Close-up view of the evolution of the normalized 
integrated intensity of the (113) peak upon pressure. Note the recovery of the pristine width and 
maximum intensity after depressurization. 



In addition to the off-centering of the Ca atoms, the study by Csanyi et al. [22[ also predicts 
a P-induced buckling of the C layers, thus leading to a lowering of the crystal symmetry. Our 
data do not enable to verify this prediction, for this would require detecting weak extra peaks, 
probably below the threshold of our diffractograms. Though, an important insight into the high- 
P structure is gained by analyzing the new features appearing above 9 GPa. These features 
must be intrinsic, for they disappear upon depressurization. Most of them are extremely broad, 
which gives evidence of an incipient amorphous phase. The only changes upon depressurization 
concern the relative intensity of some peaks. Specifically, the doublet at 2d « 9.5° (the higher 
angle peak of the doublet being the (104) reflection) becomes stronger, whilst the (003), (036) 
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and (223) peaks almost vanish. The following two phenomena may account for these changes: 
1. a P-induced texture change, which typically is irreversible; 2. a partial transformation of the 
pristine R3m phase into the competing hexagonal PQ 3 /mmc one, which would require a modest 



energy difference js]. Indeed, the diffraction patterns of the two phases are almost identical, 
but the peak intensities are different. For example, the (002) reflection of the P6 3 /mmc phase, 
which corresponds to the (003) one of the R3m phase, is very weak, in agreement with our 
observation. High resolution data on single crystals would help to verify the above scenarios. 
Finally, our results rule out the scenario of staging transition, which is known to occur in other 
graphite intercalated compounds at P ~ 0.1 GPa [231 ] . 

In conclusion, the present high pressure study shows no change of space group symmetry in 
CaCe up to 13 GPa. The experimental P-dependence of the unit cell parameters is in good 
agreement with DFT calculations. At 9 GPa, we found a jump of the isothermal compressibility 
accompanied by a large Bragg peak broadening. These two results give evidence of a P- 
induced order-disorder phase transition. Space group symmetry considerations supported by 
ab initio calculations suggest that the transition should involve the Ca sublattice. This scenario 
provides a full account of the large T c reduction and of the large increase of residual resistivity 
reported previously. Finally, the present study unveils a rather unique type of superconductor- 
superconductor transition driven by an order-disorder structural transition. 
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